arXiv:1502.05574vl [cond-mat.quant-gas] 19 Feb 2015 


POSITIVE AND NEGATIVE MASS SOEITONS 
IN SPIN-ORBIT COUPEED BOSE-EINSTEIN CONDENSATES 


V. ACHILLEOS', DJ. FRANTZESKAKIS P.G. KEVREKIDIS ^ P. SCHMELCHER 

J. STOCKHOFE^ 

'Department of Physics, University of Athens, Panepistimiopolis, Zografos, Athens 157 84, Greece 
^Department of of Mathematics and Statistics, University of Massachusetts, 

Amherst, MA 01003-9305, USA 

^ Zentrum fiir Optische Quantentechnologien, Universitat Flamburg, 

Lumper Chaussee 149, 22761 Flamburg, Germany 
"'Xhe Flamburg Centre for Ultrafast Imaging, Lumper Chaussee 149, 22761 Hamburg, Germany 
* Corresponding author E-mail: dfrantz@phys.uoa.gr 

Compiled February 20, 2015 


We present a unified description of different types of matter-wave solitons that 
can emerge in quasi one-dimensional spin-orbit coupled (SOC) Bose-Einstein conden¬ 
sates (BECs). This description relies on the reduction of the original two-component 
Gross-Pitaevskii SOC-BEC model to a single nonlinear Schrddinger equation, via a 
multiscale expansion method. This way, we find approximate bright and dark soliton 
solutions, for attractive and repulsive interatomic interactions respectively, for differ¬ 
ent regimes of the SOC interactions. Beyond this, our approach also reveals “negative 
mass” regimes, where corresponding “negative mass” bright or dark solitons can exist 
for repulsive or attractive interactions, respectively. Such a unique opportunity stems 
from the stmcture of the excitation spectmm of the SOC-BEC. Numerical results are 
found to be in excellent agreement with our analytical predictions. 

Key words: Solitons, multicomponent condensates, spin-orbit coupling. 


1. INTRODUCTION 


The recent experimental realization of spin-orbit coupling (SOC) of neutral 
atoms in Bose-Einstein condensates (BECs) |[T]-j^ and fermionic gases paved 
the way for a new era of studies in cold atom physics. Indeed, these recent develop¬ 
ments suggested the opportunity that many fundamental phenomena of condensed- 
matter or even high-energy physics can be simulated in cold atom systems. Striking 
features such as Zitterbewegung have already been observed in experiments ||^|^, 
while proposals for the study of quantum spin Hall effect |[8j, topological quantum 
phase transitions in Eermi gases Q, or even chiral topological superfluid phases 1101, 
have also been reported. It is thus clear that the perspectives are very promising: 
indeed, currently available experimental techniques allow the creation of versatile 
effective gauge potentials in ultracold atoms. Abelian or even non-Abelian, relevant 
not only to condensed-matter physics but also to the understanding of interactions 
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between elementary particles (for relevant reviews on this subject cf. Refs. 111 121). 

Although many effects related to SOC-BECs can be explored and understood 
in the linear (non-interacting) limit, there has also been much interest in the study 
of purely nonlinear effects and, in particular, in localized nonlinear excitations that 
may be supported by SOC. This way, topological structures such as skyrmions p^ , 
vortices 03. and Dirac monopoles | [T5| , as well as coherent structures in the form 
of dark 116 171 or gap 118 1, and bright |[W 201 solitons (for repulsive and attractive 
interatomic interactions, respectively) have been studied in SOC-BECs. Solitonic 
structures have already been investigated extensively in single- and multi-component 
BECs (see, e.g., the reviews | |2Tp2| ); nevertheless, the above mentioned recent stud¬ 
ies have shown that the presence of SOC enriches significantly the possibilities re¬ 
garding the structural form of solitons, as well as their stability and dynamical proper¬ 
ties. As an example, we note the recent work |23 1, where beating dark-dark solitons, 
that occupy both energy bands of the spectrum of a SOC-BEC and may perform Zit- 
terbewegung oscillations, were predicted to occur. Interestingly, the solitonic struc¬ 
tures studied in Ref. p^ , are embedded families of bright, twisted or higher excited 
solitons inside a dark soliton which, as well, are only possible due to the structure 
of the energy spectrum of SOC-BECs. Moreover, phenomena stemming from the 
interplay of SOC with additional manipulation capabilities available in BECs such 
as optical lattices or Raman coupling fields are presenfly under infense investigation. 
Canonical examples fhereof can be seen in fhe manifesfafion of dynamical insfabil- 
ifies of excifafions in fhe vicinify of a band gap |24| of optical laffices, or in fhe 
demonsfrafion of a Dicke-fype quanfum phase fransifion, by changing fhe Raman 
coupling sfrengfh befween a spin-polarized and a spin-balanced phase, in analogy fo 
quanfum optics p5] |. 

If is fhe purpose of fhis work fo provide a unified descripfion of differenf fypes 
of solitons fhaf may occur in eifher ahracfive or repulsive quasi one-dimensional (ID) 
SOC-BECs. Our mefhodology relies on fhe sfudy of solifon solufions of an effecfive 
scalar nonlinear Schrbdinger (NES) equafion: fhis can be derived from fhe sysfem 
of fhe fwo coupled Gross-Pifaevskii equafions describing fhe SOC-BEC via a mul- 
fiscale expansion mefhod. We firsf find eifher dark or brighf solifon solufions, for 
repulsive or ahracfive inferacfions, respecfively, of amplifudes and widfhs depending 
on fhe SOC paramefers. These solifons may be formed in eifher of fhe fwo regions of 
fhe excifafion energy specfrum, where fhe lowesf energy band has eifher a single- or 
a double-well sfrucfure; so-called “sfriped” solifons, in fhe region where fhe double¬ 
well of fhe specfrum is symmefric, are shown fo be possible, foo. The above men¬ 
tioned fypes of solifonic sfrucfures will be called hereaffer “positive mass solifons”, 
fo distinguish fhem from fhe possibilify of “negafive mass” ones infroduced in fhe 
following. 

Indeed, inferesfingly enough, our analyfical approach reveals fhe exisfence of 
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“negative mass” parametric regions, meaning that the sign of the kinetic energy term 
(the dispersion coefficient in the effective NLS model) becomes inverted. In such 
negative mass regions, we show that bright or dark solitons can be formed in SOC- 
BEC with repulsive or attractive interactions, respectively. Pertinent solitonic struc¬ 
tures will be called hereafter “negative-mass solitons”. Notice that such solitons can¬ 
not exist in usual single-component BECs described by a scalar NES model: there, 
bright (dark) matter-wave solitons are formed only for BEC with attractive (repul¬ 
sive) interactions, respectively p^ . Instead, bright solitonic structures, for instance, 
can only exist as gap solitons in repulsive BECs loaded in an optical lattice potential, 
where the effective mass also changes sign (becomes negative) at the edge of the first 
Brillouin zone [271; for a theoretical proposal of the extension of such states in SOC- 
BECs, see, e.g., 1181. Thus, our results suggest an alternative method of “dispersion 
management” of matter-waves based on the SOC effect (as previous proposals relied 
on the use of optical lattices |[28|). 


2. THE MODEL AND ITS ANALYTICAL CONSIDERATION 

We consider the case of a quasi-lD SOC-BEC, confined in a frap wifh longifu- 
dinal and fransverse frequencies, Ux and a;_L, such fhaf <C tuj,. In fhe framework 
of mean-field fheory, fhis sysfem is described by fhe energy funcfional : 

= + ^ (5fll|V’l|^ + fl'22|V’2|‘^ + 2fl'12|V’l|^IV'2p) , (1) 

where 'J' = (V’i,'02)^> and fhe condensate wavefunctions V’ 1,2 are the two pseudo¬ 
spin components of the BEC. Eurthermore, the single particle Hamiltonian Tfo in 
Eq.Q reads: 


'Ho = ^{px + kL&z)‘^ + Vtr{x) + nax + 6az, (2) 

2m 

where px = —ihdx is the momentum operator in the longitudinal direction, m is the 
atomic mass, and ax^z are the Pauli matrices. The SOC terms are characterized by 
the wavenumber ki of the Raman laser which couples the two components, and the 
strength of the coupling O. Moreover, S is an energy shift due to the detuning from 
Raman resonance. The external trapping potential Vtr{x), is assumed to be of the 
usual parabolic form, Efr = (l/2)ma;^x^. Einally, the effective ID coupling con¬ 
stants Qij are given by gij = 2 hL 0 ±aij, where Uij are the s-wave scattering lengths; 
below, we will present results for both attractive (oij < 0) and repulsive {aij > 0) 
interatomic interactions. 

Using Eq. we can obtain the following dimensionless equations of motion: 
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= ^-^5^-iA:L(9^ + Vtr + <7(|V’i|^ + /3|V'2p)+<5^ V’i + f^V'2, (3) 

= {^]^dl + ikLdx + Vtr + cr{f3\'ilJi\^+ \'iIj 2?)-^'ip2 + ^i^u (4) 

where energy, length, time and densities are measured in units of a_\_ (which 
is equal to ^Jh/muj±), cuj^, and an, respectively, and we have also used the trans¬ 
formations ^ a±kL, ^ ^ Q/{huj±) and 6 —)• d/{huj±). Notice that we have also 

assumed that the intraspecies scattering lengths are identical, an = a 22 ; this way, 

parameter a, defined as cr = sgn(an) = sgn(a 22 ) = ±1, corresponds to repulsive 
and attractive intra-species interactions, respectively, while /3 = ai 2 /|an|. Finally, 
the trapping potential in Eqs. §-S is now given by Vtr{x) = {\/2)uj‘^^x^, where 
Wtr = Waj/tux ^1- To a first approximation, when sufficiently weak, the trapping 
potential can be neglected. Thus, below we will only consider the homogeneous 
case, i.e., Vtr = 0 (nevertheless, we have checked that the results do not change qual¬ 
itatively upon including the parabolic trap - see also Refs. (nKig). 

Below we will use the following parameter values: k^ = 8 for the Raman 
wavenumber, fl G [0,120] for the normalized Raman coupling strength, which are 
relevant to the experiments of Ref. Q. Finally, for our analytical considerations, /3 
will be treated as an arbitrary (positive) parameter. In fact, for the hyperfine states of 
®^Rb commonly used in today’s SOC-BEC experiments this parameter takes values 
/? 1, due to the fact that the inter- and intra-species scattering lengths take ap¬ 

proximately the same values. Eor this reason, in the simulations that we will present 
below, this parameter is fixed to the value /3 = 1; nevertheless, we have confirmed that 
physically relevant deviations from this value do not alter qualitatively the results. 

We now use a multiscale perturbation method to reduce the system of Eqs. Q- 
Q to a scalar NES equation. Such a reduction will allow us to derive approximate 
analytical bright and dark soliton solutions, depending on the type of the interactions. 
To proceed, we introduce the ansatz: 


oo oo . ^ 

n=l n=l ^ ^ / 

where the vectors u„ = [[/„, Vn]'^(j)n are composed by the coefficients Un and Vn and 
the unknown field envelopes (jin = 4’n{T,X). The latter are assumed to be functions 
of the slow variables T = and X = e{x — vt), where v is the velocity (to be 
determined). Additionally, k is the momentum, p. = u} + e^ujQ the chemical potential; 
here, u is the energy in the linear limit, e^wo is a small deviation about this energy 
(e <C 1), and ojq/uj = 0(1). Introducing Eq. Q into Eqs. ([^-Q, we derive the 
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Fig. 1 - (Color online) Sketches of the linear energy spectrum u)(k) for Region 1 and Region 2, 
denoted by R 1 and R 2; cf. left and right panels, respectively. Top (bottom) panels correspond to 
(5 = 0 (<5 7^ 0), and dashed lines depict the group velocity dispersion Lj''{k) of the lower band. 

following set of equations at orders 0{e), O(e^) and 0{e^), respectively: 


Wui = 0, (6) 

Wu2 = iWo^xui, (7) 

Wu3 = iWo5xU2 - (^dT + ]^d\-A + oj^ ui, (8) 

where Wq = (W — loI)' (primes denote hereafter differentiation with respect to k), 
while matrices W and A are: 


ui — k'^/2 — kkL — S —Q 

—n uj — k'^/2 +kkL +6 ’ 

r7(|[/i|2+/3|yi|2) 0 1 

0 a(/3|i7i|2 + |yi|2) _ • 


(9) 

( 10 ) 


At 0(e), the solvability condition detW = 0 yields the linear excitation energy 
spectrum 


UJ = U±{k) = h‘^± ^J{kkL + <5)2 + (72, (11) 

which consists of two different bands, the upper and lower one, as shown in Fig.[T] 
The lower band of the spectrum uj- has a single local minimum for Q./k'j^> 1, 
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and two local minima iox Q./k\ < 1; respective regions of the parameter space will 
be called hereafter “Region 1” and “Region 2” respectively (cf. left and right panels 
of Fig. [T]l. Notice that for (5 = 0 the spectrum is symmetric with respect to /c = 0, 
while for (5 7 ^ 0 this symmetry is broken (cf. top and bottom panels of Fig. [T]). We 
should also mention, at this point, that below we will focus on nonlinear states, in 
the form of matter-wave solitons, corresponding to the lower energy band (ground 
state). However, as shown in Ref. p^ , solitons may also exist in an excited state 
occupying simultaneously both energy bands | [^ . Finally, at the same order 0{e), 
the compatibility condition of Eq. Q yields: 


ui = ( ^ 

where the parameter Q is given by: 


( 12 ) 


Q = Q{uj,k) 



kki-S). 


(13) 


Notice that the above parameter sets the left and right eigenvectors of W at 
eigenvalue 0, given by L = [1, Q] and R = [1,(5]^, respectively. 

Next, at the order 0{e^), the compatibility condition of Eq. ([^, LWqR = 0, 
fixes V at the group velocity, v = oj'{k), yielding the explicit condition 


0^-1 

^ = (14) 

Q‘^ + 1 

This means that solutions for (j) in the lower energy band may be either static 
for k = /cjnin (choosing a minimum kmin of the dispersion where the group velocity 
vanishes), or moving for k ^ /cmin- Note that, in the following, both static and moving 
soliton solutions will be presented. At the same order, we also obtain the form of the 
solution for U 2 : 


^dxMX,T). (15) 

Einally, the compatibility condition for Eq. ([^ at 0{e^), together with Eqs. ( [T^ 
and ( [T5] ), yields the following scalar NES equation for the unknown field (f)i\ 

idT(l) =d\(f) +—ojQcj). (16) 

where fhe subscripf has been omiffed for convenience, while paramefer v is given by: 



v{k) = 


g4 + 2/3Q2 + i 


1 + Q2 


(17) 
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and it is always positive, since parameter Q [cf. Eq. ( [T3] )] is real (note that here we 
have assumed that (3 > 0). The above NLS model plays a central role in our analysis, 
as it will be used below for the derivation of various soliton families. 


3. SOLITON SOLUTIONS 


The sign of the nonlinearity and dispersion coefficients of the above effective 
NLS model is crucial in determining the types of solitons that can be supported: for 
u''a > 0 the solitons are dark, while for a;"f7 < 0 the solitons are bright. Obviously, 
the sign of the nonlinearity coefficient depends only on the type of interatomic inter¬ 
actions: cj = ±1 for repulsive and attractive interactions, respectively. 

On the other hand, let us take a closer look at the dispersion coefficient oj" = 
d'^Lo/dk'^, which can a lso be viewed as the inverse of the effective mass, i.e., uo" oc 


(see, e.g.. Refs. |27 ^). The dependence of this coefficient on momentum 
k is depicted by the dashed lines in Fig. [T] It is observed that in the lower energy 
band, in both Regions 1 and 2, and for k = femin (where group velocity is zero, as 
mentioned above), one has always uj" > 0. This indicates the existence of “positive 
mass” stationary dark (bright) solitons in SOC-BECs with repulsive (attractive) in¬ 
teractions. Positive mass moving solitons can exist as well, featuring a finite group 
velocity that can be found for k ^ /cmin- Notice that dark solitons exist inside the 
linear band (wq > 0), while the bright solitons are found in the infinite gap below 
the lower energy band (wq < 0); this can be readily seen by the stationary form of 
Eq. ( [T^ , corresponding to dTfj) = 0 

On the other hand, in Region 2, there exist intervals of momenta (depicted by 
the ellipses in the right panels of Fig. where oj” < 0. Thus, also “negative mass” 
moving dark (bright) solitons can exist in SOC-BECs with attractive (repulsive) in¬ 
teractions. 

We now proceed by presenting different types of exact soliton solutions of 
Eq. ( [T^ and corresponding approximate soliton solutions of the original system of 
Eqs. Q-©- We will also study the solitons’ properties by means of direct numerical 
simulations. 


3.1. POSITIVE MASS SOLITONS 


We Start with the case of positive mass solitons, which may be either stationary 
(exhibiting a momentum k = fcmin at the minima of the energy spectrum where the 
group velocity v vanishes) or moving (for momenta k ^ /cmin)^ as mentioned above. 

For repulsive interactions, a > 0, we may write the dark soliton solution of 
Eq. ( [T^ , with Wo > 0, in the form |21 261: 
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0d = (cos discnhzd + i sin6*) 


(18) 


where Zd = -s/uJoJuf cos9[X — Xo(T)]; here, Xq{T) is the soliton center, while the 
so-called soliton phase angle 9 controls the soliton amplitude, -^JoJ^^cosO, and soli¬ 
ton velocity in the slow time scale T through the equation dX^/dT = s/ ojq/ uj" smO. 
Notice that the above soliton solution is characterized by two free parameters, one 
for the background, uiq, and one for the soliton, 0. 

For attractive interactions, u > 0, a bright soliton solution of Eq. ( [T^ , with 
ojQ < 0, may be written as: 


0b = r]sec]iZhexp{iKX) 


(19) 


where Zb = u/oj''[X — Xo(T)]; here, rj is the soliton amplitude which is con¬ 
nected with the parameter loq through the relation coq = ju") j2, Xq(T) 

is the soliton center, while the wavenumber n is connected with the soliton velocity 
through dXo/dT = uj''k. As in the previous case, the above soliton solution has two 
free parameters, ujq and k. 

Having found the solutions of Eq. we may write the approximate sta¬ 
tionary soliton solutions for the two components of the SOC-BEC, in the original 
coordinates as follows: 



for the dark and bright soli tons respectively; recall that Q depends on k [or w, as per 
the dispersion relation ([TT])]. 

Another family of soliton solutions can also be found as follows. In Region 2, 
and for (5 = 0 (see top right panel of Eig. [T]), there exist two degenerate minima at 
specific values of momenta, say kc and —kc, for which two different nonlinear solu¬ 
tions can be obtained, as per the asymptotic analysis presented above. However, in 
the linear limit, a superposition of these states is also a solution with the same energy. 
Using a continuation argument, it turns out to be possible to find a nonlinear solution 
as well. In the case <5 = 0 and /3 = 1, we may employ the symmetry of Eqs. Q-Q, 
namely 0i = —02 (bar denotes complex conjugation), and find approximate solu¬ 
tions of the original equations in the form: 
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Fig. 2 - (Color online) Positive mass dark solitons for ct = -|-1. In the top panels shown are the real 
parts of '01,2 1 while the bottom panels depict the densities of the two components. The left column 
corresponds to Region 1 with = 1.56 and the middle panel to Region 2 with = 0.78 both 
for (5 = 20. The right panels depict the “stripe” solitons for = 0.78 and <5 = 0. Dashed (hlack) 
lines in the bottom panels depict corresponding analytical results. 


f ^ ^ rrh f Q+^os{kcx) + iq-sm{kcx) \ 

V V ’2 y ^ V —Q+cos{kcx) + iq-sm{kcx) y’ 

where q± = + Q{±kc) and (7 is a free parameter. 

Examples of all three different families of dark and bright soliton solutions are 
shown in Figs. and respectively. These states were obtained numerically, by 
solving (via a fixed point algorithm) the stationary form of Eqs. ([^-(Q obtained by 
using = $(x) exp(—f/if); then, for $(a:) the stationary form of the solitons 

of Eqs. ( [20| ), ( [2T] ) and ( [22] ) was used as an initial guess. Notice that for obtaining 
these solitons, we have used the momentum k = /cmin^ corresponding to zero group 
velocity. In other words, the solitons shown in Figs. [^ and [^ are stationary ones; for 
arbitrary momenta, positive mass moving solitons also exist, but we will not discuss 
them further (they were found to exist and be stable in parameter regions where their 
stationary counterparts are stable 

In the left panel of Fig. [^ (Fig. [^, a dark (bright) soliton in Region 1 is shown 
for ^/k\ = 1.56 and 6 = 20. In the middle panel of these figures, we show respec¬ 
tive solitons in Region 2, for fi/k^ = 0.78 and S = 20. Note that the number of 
atoms of the two components is generally different, with the asymmetry arising due 
to the presence of the constant factor Q [cf. Eqs.(|20|)-(|2T])]. Finally, in the right 
panel of Fig. [^ (Fig. [^, an example of a “stripe” dark (bright) soliton is shown for 
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Fig. 3 - (Color online) Positive mass bright solitons for cr = — 1. In the top panels shown are the real 
parts of '01,25 while the bottom panels depict the densities of the two components. The left column 
corresponds to Region 1 with ^/k\ = 1.56 and the middle panel to Region 2 with = 0.78 both 
for 5 = 20. The right panels depict the “stripe” solitons for ^/k\ = 0.78 and <5 = 0. Dashed (hlack) 
lines in the bottom panels depict corresponding analytical results. 


VLjk\ = 0.78 and J = 0. Note that in this case the two components have the same 
amplitudes and number of atoms, which is a consequence of the symmetry of these 
states, associated with the fact that yii = — 1 /’ 2 - In both figures (Figs. and |^, we 
also plot the analytical result of Eqs. (|20|)-( 221, depicted by a dashed (black) line, 
illustrating the very good agreement between our approximate analytical result and 
the numerically obtained solutions. 


3.2. NEGATIVE MASS SOLITONS 

We now proceed with the case of negative mass solitons, which can be found in 
specific domains of Region 2, where a;" < 0 (cf. ellipses in fhe righf panels of Fig. [^. 
In such a case, if is possible fo obfain brighf solifon solufions for a repulsive SOC- 
BEC, and dark solifon solufions for an attractive BEC. This is a special fealure of fhe 
spin-orbif coupling and of fhe resulting sfrucfure of fhe linear energy specfrum. In 
facl, fhis is a fealure fypically absent from homogeneous (especially single compo- 
nenl) syslems, allhough such a possibilily can be generaled in inhomogeneous BECs, 
e.g. in fhe presence of an oplical lallice p7| . Negative mass brighf and dark solifon 
solutions have fhe same functional form as positive mass ones [Eqs. (|20|)-(|2T])], bul 
now for cr = +1 and cr = — 1, respectively. 

In order lo illusfrale fhe exislence and dynamics of such sfales, we numerically 
inlegrale Eqs. (|^-(0 wilh an inilial condition corresponding lo a brighf solifon in 
fhe repulsive case and lo a dark solifon for fhe alfraclive case. The evolution of a 
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Fig. 4 - Top panel: contour plot showing the evolution of the total density I'i/'iP -I- |'!/' 2|^5 for an initial 
bright soliton in the case of repulsive interactions (a — Bottom panel: same as above, but now 
for a dark soliton in the case of attractive interactions (a = —1). Parameter values used in both cases 
are 0./k\ = 0.49, 5 = A,k = —2.3 and = 0.1. 


bright soliton propagating in a SOC-BEC with repulsive interactions, is shown in the 
top panel of Fig. with parameter values Q./k\ = 0.49, <5 = 4, e^tuo = 0.1 and a 
momentum k = —2.3. On the other hand, the evolution of a dark soliton, but for 
an attractive SOC-BEC, is shown in the bottom panel of the same figure, with the 
same parameter values as above. In both cases the solitons are found to be robust, 
and propagate without any distortion or significant emission of radiation; in fact, we 
have confirmed their robustness for times up to at least t = 400 in dimensionless 
units. Notice that in the numerical simulations we have also checked that the plane 
wave background of a negative mass dark soliton (cr = — 1) is not subject to the 
modulational instability: in the simulations, we used a perturbation of wavenumber 
K outside the negative-mass region (e.g., with K = /cmin) and found that the soliton 
background remains modulationally stable. 

The negative mass solitons shown in Fig. feature a relatively large velocity: 
for a typical parameter value of the transverse frequency u!± ^ 27r x 400 Hz for 
Rubidium atoms, the solitons propagate about 65 microns in 0.4s (i.e., the velocity is 
|u| = 162 fim/s or |u| = 1.2 in dimensionless units). Nevertheless, slower (and even 
stationary) solitons can also be obtained upon properly choosing different momenta 
values. A pertinent example is shown in the top panel of Fig.[^ there, a bright soliton 
of momentum k = —0.9, corresponding to a group velocity v = —0.05, is depicted. 

For the above case of slower solitons, and to better illustrate the fact that nega¬ 
tive mass solitons exist only due to the presence of spin-orbit coupling, we perform 
the following numerical experiment (cf. Fig.[^. At time t = 90, we turn off the SOC 
terms upon setting Q,kL and 6 equal to zero. Then, as it is expected for a conden- 
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Fig. 5 - (Color online) Same as in Fig.|^but for k = —0.9. The vertical dashed lines indicates the 
instant when the spin-orbit coupling is switched off and the insets show each component’s density at 

f = 10 (left) and t = 110 (right). 


sate with repulsive interactions, the bright soliton is not supported and, as a result, it 
begins to spread. The insets show the density profile of the two components 
when spin-orbit coupling is present (left) and when it is turned off (right). 

We have also studied the case of a negative mass dark soliton, propagating in 
an attractive BEC. A pertinent example is shown in the bottom panel of Fig. for the 
same parameter values (the dark soliton also has the same group velocity as its bright 
counterpart). After the initial undisturbed propagation of the soliton, we again switch 
off the spin-orbit coupling, at f = 90. Then, as it is natural for the attractive case, 
the dark soliton can no longer be supported and is destroyed (its associated density 
dip spreads), while the background on which it was supported becomes subject to a 
modulational instability, gradually producing an array of bright solitary waves. 

Note that without switching off the SOC terms, we have confirmed the robust 
evolution of these “slow” negative mass bright and dark solitons for the repulsive and 
attractive case, respectively, for up to f = 400 in dimensionless units. 


4. CONCLUSIONS 

Concluding, in this work we have provided a unified description of different 
types of solitons that may occur in quasi-ID SOC-BECs for both repulsive and at¬ 
tractive interactions. Our approach relies on a multiscale expansion method, which 
was used to reduce the original system of two coupled Gross-Pitaevskii equations 
into an effective scalar nonlinear Schrddinger (NFS) equation. 

Investigating the dispersion coefficient of the effective NFS model, we identi¬ 
fied regions where it can become either positive or negative. This important effect. 
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which arises only due to the presence of spin-orbit coupling, is reminiscent of what 
happens in the case of single-component BECs trapped in optical lattices: in the lat¬ 
ter setting, due to the possibility of change of sign of the dispersion (effective mass 
in the framework of the relevant Gross-Pitaevskii model), bright (gap) solitons can 
be formed in BEC with repulsive interactions p7| . This suggests the use of SOC 
for dispersion management of matter-waves, as in the case of BECs loaded in optical 
lattices | |28l . 

Employing the above finding, we presented both positive and negative mass 
solitons that can be formed for positive and negative dispersion coefficient, respec¬ 
tively. In particular, positive mass dark (bright) solitons were found for repulsive 
(attractive) interactions, in either of the two regions of the energy spectrum, where 
the lowest energy band has either a single- or a double-well structure. Additionally, 
it was shown that “striped” solitons, in the region where the double-well of the spec¬ 
trum is symmetric, are possible, too. These positive mass solitons may exist in either 
stationary or moving form. 

On the other hand, we have shown that negative mass bright (dark) solitons can 
be formed in SOC-BEC with repulsive (attractive) interactions; these solitons can 
also exist either as stationary or moving objects. Using direct numerical simulations, 
we found that such negative mass solitons are quite robust and can only exist in the 
presence of spin-orbit coupling. Indeed, letting the solitons evolve up to a certain 
time and then switching off the SOC terms, we found that they are not supported by 
the system without these terms. 

Our methodology and findings suggest many interesting future research direc¬ 
tions. Eirst, it would be interesting to further investigate and propose mechanisms of 
SOC-mediated dispersion management of matter-waves. Additionally, it would be 
relevant to generalize our results for multi-component, F = 1 spinor SOC-BEC 1311 
and study solitons in such systems. Another interesting theme deserving attention 
is the analysis of modulational instability of plane waves in the original system of 
Eqs. such an analysis, apart from being interesting by itself, would be par¬ 

ticularly relevant to the stability of dark solitons in this setting. Einally, it would 
be particularly interesting to extend the above considerations to higher-dimensional 
settings and explore the dynamics, e.g., of vortices or vortex rings in the presence of 
such effective dispersions. 

Acknowledgements. The work of V.A. and D.J.E. was partially supported by 
the Special Account for Research Grants of the University of Athens. V.A., D.J.E. 
and RG.K. appreciate warm hospitality at the University of Hamburg, where part of 
this work was carried out. RG.K. acknowledges support from the National Science 
Eoundation under grants CMMI-1000337, DMS-1312856, from ER7-Reople under 
grant IRSES-605096, from the Binational (US-Israel) Science Eoundation through 
grant 2010239 and from the US-AEOSR under grant EA9550-12-10332. J.S. ac- 




14 


V. Achilleos et al. 


knowledges support from the Studienstiftung des deutschen Volkes. 


REFERENCES 

1. Y.-J. Lin, K. Jimenez-Garcia, and I. B. Spielman, Nature, 471, 83 (2011). 

2. V. Galitski, and I. B. Spielman, Nature, 494, 54 (2013). 

3. J.-Y. Zhang el al, Phys. Rev. Lett. 109, 115301 (2012). 

4. P. Wang et al, Phys. Rev. Lett. 109, 095301 (2012). 

5. L. W. Cheuk et al, Phys. Rev. Lett. 109, 095302 (2012). 

6. C. Qu et al, Phys. Rev. A 88, 021604(R) (2013). 

7. L. J. LeBlanc et al. New J. Phys. 15, 073011 (2013). 

8. C. J. Kennedy et al, Phys. Rev. Lett. Ill, 225301 (2013). 

9. M. Gong, S. Tewari, and C. Zhang, Phys. Rev. Lett. 107, 195303 (2011). 

10. X.-J. Liu, K. T. Law, and T. K. Ng, Phys. Rev. Lett. 112, 086401 (2014). 

11. J. Dalihard, F. Gerbier, G. Juzeliunas, and P. Ohberg, Rev. Mod. Phys. 83, 1523 (2011). 

12. N. Goldman, G. Juzeliunas, P. Ohberg, 1. B. Spielman, Rep. Progr. Phys., 77, 126401 (2014). 

13. T. Kawakami, T. Mizushima, M. Nitta, and K. Machida, Phys. Rev. Lett. 109, 015301 (2012); 

C.-F. Liu and W. M. Liu, Phys. Rev. A 86, 033602 (2012); 

X. Zhou, Y. Li, Z. Cai, and C. Wu, J. Phys. B: At. Mol. Opt. Phys. 46, 134001 (2013). 

14. X.-Q. Xu and J. H. Han, Phys. Rev. Lett. 107, 200401 (2011); 

J. Radic T. A. Sedrakyan, 1. B. Spielman, and V. Galitski, Phys. Rev. A 84, 063604 (2011); 

B. Ramachandhran et al, Phys. Rev. A 85, 023606 (2012); 

A.L. Fetter, Phys. Rev. A 89, 023629 (2014). 

15. G. J. Conduit, Phys. Rev. A 86, 021605(R) (2012). 

16. O. Fialko, J. Brand, and U. Zulicke, Phys. Rev. A 85, 051605(R) (2012). 

17. V. Achilleos, J. Stockhofe, P. G. Kevrekidis, D. J. Frantzeskakis, and P. Schmelcher, Europhys. 

Lett. 103, 20002 (2013). 

18. Y. V. Kartashov, V. V. Konotop and F. Kh. Abdullaev, Phys. Rev. Lett. Ill, 060402 (2013); 

V. E. Lobanov, Y. V. Kartashov, and V. V. Konotop, Phys. Rev. Lett. 112, 180403 (2014). 

19. V. Achilleos, D. J. Frantzeskakis, P. G. Kevrekidis, and D. E. Pelinovsky, Phys. Rev. Lett. 110, 

264101 (2013). 

20. Y. Xu, Y. Zhang, and B. Wu, Phys. Rev. A 87, 013614 (2013). 

21. D. J. Frantzeskakis, J. Phys. A 43, 213001 (2010). 

22. Y. Kawaguchi and M. Ueda, Phys. Rep. 520, 253 (2012). 

23. V. Achilleos, D. J. Frantzeskakis, and P. G. Kevrekidis, Phys. Rev. A 89, 033636 (2014). 

24. C. Hamner et al, arXiv: 1405.4048 (2014). 

25. C. Hamner et al, Nat. Commun. 5, 4023 (2014). 

26. P. G. Kevrekidis, D. J. Frantzeskakis, and R. Carretero-Gonzalez, Emergent Nonlinear Phenom¬ 

ena in Bose-Einstein Condensates (Springer-Verlag, Berlin, 2008); R. Carretero-Gonzalez, D. J. 
Frantzeskakis, and P. G. Kevrekidis, Nonlinearity 21, R139 (2008). 

27. B. Eiermann et al, Phys. Rev. Lett. 92, 230401 (2004). 

28. B. Eiermann et al, Phys. Rev. Lett. 91, 060402 (2003). 

29. T. L. Ho and S. Zhang, Phys. Rev. Lett. 107, 150403 (2011); 

S. Sinha, R. Nath, and L. Santos, Phys. Rev. Lett. 107, 270401 (2011); 

Y. Li, L. P. Pitaevskii, and S. Stringari, Phys. Rev. Lett. 108, 225301 (2012). 



Positive and negative mass solitons in spin-orbit coupled Bose-Einstein condensates 


15 


30. The excitation of the upper band was demonstrated in the experiments of Refs. (6|7| , where Zitter- 

bewegung oscillations were observed. Similarly, solitons occupying both upper and lower bands 
may, under certain conditions, exhibit Zitterbewegung oscillations as well |23| . 

31. Z. Lan and P. Ohberg, Phys. Rev. A 89, 023630 (2014). 



